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Abstract 

The International Muon Ionization Cooling Experiment (MICE) is devoted to a 
study of a muon cooling channel capable of giving the desired performance for a Neu- 
trino Factory. One of the goals is achieving an absolute accuracy of measurements of 
emittance reduction as high as ±0.1%. The paper describes results of a Monte Carlo 
study on allowed density variations of liquid hydrogen corresponding to the desired 
accuracy of the measurements. 

1 Introduction 

An experiment which allows the investigation of the performance of a muon cooling chan- 
nel at different conditions of interest has been designed by MICE collaboration Q. In par- 
ticular, high precision measurements of emittance reduction are planned to be performed 
with an accuracy of ±0.1%. Various factors should be taken into account to achieve the 
desired accuracy. One of the essential parameters is material density of the liquid hydrogen 
in the absorber. The density is influenced by temperature and pressure of the hydrogen kept 
within pre-determined tolerances. The Monte Carlo study with the MARS 14 code was 
undertaken to calculate the dependence of the emittance reduction on hydrogen density at 
realistic conditions. The results can be used to determine the above-mentioned tolerances 
corresponding to the desired accuracy of measurements. 

2 Basic Formulae 

In a four-dimensional phase space with coordinates x, p x , y, p y , the normalized emittance 
of a muon beam, e n , can be calculated according to the expression El 
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where Z is a 4 x 4 correlation matrix, 
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and a is x — x, b is (p x — p x )/m jJ c, x and /? x are muon coordinate and momentum, re- 
spectively, along x— axis. The entries c and rf are for y— axis and analogous to a and 
respectively. In the expressions the bar over the symbols means statistical averaging over 
an ensemble of simulated muon trajectories. 

To generate a proper muon beam distribution in the system, symmetry considerations 
are taken into account H. Firstly, Gaussian distributions are modeled for x,p x ,y,p y in the 
geometrical center of the absorber, where magnetic field flips, using information on the 
P-function distribution in the channel Q. Secondly, backward muon transport is modeled 
in the magnetic field and without the material (hydrogen and aluminum) to get a reflected 
incoming muon distribution in front of the absorber. And thirdly, regular muon transport 
through the absorber is modeled. 

For a given ensemble of trajectories, the matrix Z is calculated both in front of and 
behind the liquid hydrogen absorber. After that the emittance reduction, Ae n , is calculated 
according to expression (0Q). 



Realistic geometry of the 35-cm liquid hydrogen absorber along with 'inflected' window 
design |Hl|5l was implemented in the model (see Fig.[T]). The absorber is inside a solenoid 
which is taken into account via its magnetic field only. Thickness of the aluminum windows 
of the absorber is variable with the thinnest parts being on the z-axis. 

Realistic three-dimensional distributions of magnetic field over the cooling channel 
were taken from Ref. B. The longitudinal and radial solenoidal field distributions as im- 
plemented in the MARS model are shown in Fig. [2j 



The Monte Carlo calculations were performed for 200-MeV/c incident muons at various 
densities of the liquid hydrogen. Results of the calculations are shown in Fig.[3J Statistical 
uncertainty (la) of the calculations was about 0.02% which is less than the linear size of 
the symbols used in the Figure. 



3 Geometry Model 



4 Calculated Results 
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Figure 1 : A cross section of the MARS model of a MICE liquid hydrogen absorber along 
with 100 sampled muon tracks. 
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Figure 2: The longitudinal (left) and radial (right) distributions of magnetic field around 
the liquid hydrogen absorber. The arrows indicate the field direction only, not magnitude. 
The field flips at the geometrical center of the absorber. 
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Figure 3: Calculated emittance reduction vs liquid hydrogen density in the central absorber 
of the cooling channel for 200-MeV/c muons Q. Here 100% corresponds to 0.0708 g/cm 3 . 

One can see that a variation as high as 2% in hydrogen density gives rise to a 0.1% 
variation in Az n . Therefore, to ensure the desired accuracy of emittance measurements, 
tolerance for the density of the liquid hydrogen should not exceed 2%. 

5 Conclusions 

Monte Carlo calculations were performed on emittance reduction of a muon beam vs hy- 
drogen density for a MICE liquid hydrogen absorber within realistic absorber geometry and 
with detailed three-dimensional distribution of magnetic field. It was shown that, within 
the range of interest, the density dependence is clearly linear. To ensure the accuracy of 
emittance measurements as high as 0.1%, tolerance for the hydrogen density should be less 
than 2%. 
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